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EFFECT  OF  FINITE  CURRENT  CHANNEL  WIDTH 
ON  THE  COLLISIONAL  ION  CYCLOTRON  INSTABILITY 


I.  INTRODUCTION 

It  is  well  known  that  an  equilibrium  magnetic  field-aligned  current 


can  result 

in 

the 

growth  of  obliquely  propagating  electrostatic 

ion- 

cyclotron 

waves 

in 

collisionless  [Drummond 

and 

Rosenbluth,  1962] 

and 

eollisional 

[Mine*, 

1972;  Chaturvedi  and 

Kaw, 

1975]  plasmas. 

This 

instability  has  extensively  been  investigated  in  the  laboratory  machines 
[D'Angelo  and  Motley,  1962;  Cartier  et  al.,  1985],  and  is  believed  to  have 
been  observed  in  the  auroral  ionosphere  [Kelley  et  al.,  1975;  Yau  et  al., 

1  983;  Fejer  et  al.,  1984].  The  linear  and  nonlinear  theory  of  this 
instability  applied  to  the  auroral  ionosphere  is  also  well  developed 
[Kindel  and  Kennel,  1971;  Chaturvedi,  1976;  Satyanarayana  et  al.,  1985]. 

-  The  field-aligned  currents  in  the  auroral'  ionosphere  are  an  integral  part 
of  the  vionosphere-magnetosphere  coupling  system  and  it  is  fairly  well 
established  now  that  these  currents  can  be  highly  nonuniform  in  structure. 
[Burke  et  al.,  1983;  Bythrow  et  al.,  1984],  Often,  the  current  system  is 
in  the  form  of  sheets  which  have  variable  thicknesses.  Earlier  efforts  to 
study  the  oscillation  modes  of  the  finite  current  sheets  included  those  by 
Elliott  [1975]  and  Dungey  and  Strangeway  [1976].  1  An  attempt  to  study  the 
effect  of  finite  width  of  a  current  sheet  on  the  current  driven 
collisionless  ion  acoustic  instability  was  made  by  Hwang  et  al.  [1983Pwith 
applications  to  the  auroral  situation.  They  found  that  the  finite 
thickness  of  a  current  sheet  results  in  partially  stabilizing  the  system 
and  contributes  to  the  coherence  of  the  excited  waves.  Bakshi  et  al. 
[1983]  have  studied  the  problem  of  finite  width  currents  on  the 

collisionles3  current  driven  ion  cyclotron  instability  and  found  that  for 
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sufficiently  narrow  current  channels  the  mode  is  stabilized  (for  < 
few  p  where  is  the  width  of  the  current  channel  and  p is  the  mean  ion 
Larmor  radius).  In  this  paper  we" investigate  the  effects  of  a  finite  width 
current  channel  on  the  collisional  ion  cyclotron  (CICI)  instability.  We 
find  that  the  instability  can  also  be  stabilized  for  sufficiently  narrow 
current  channels. 

The  organization  of  the  paper  is  as  follows.  In  the  next  section  we 
present  the  basic  assumptions  and  the  derivation  of  the  mode  equation.  In 
Section  III  we  present  an  analysis  of  the  mode  equation,  both  analytical 
and  numerical.  In  the  final  section  we  summarize  the  results  and  discuss 
applications  of  the  theory  to  laboratory  and  space  plasmas. 

II.  DERIVATION  OF  THE  MODE  EQUATION 

The  geometry  and  plasma  configuration  used  in  the  analysis  is  shown  in 
Fig.  1.  The  ambient  magnetic  field  is  uniform  and  in  the  z-direction 
(.B  »  3^  ez],  while  equilibrium  current  is  non-uniform  in  the  x-direction 

a 

and  is  directed  in  the  z-direction  [j  =  Jn(x)  e  }.  The  current  is  assumed 

—  u  z 

to  be  carried  by  the  electrons.  The  density  and  temperature  are  taken  to 

be  homogeneous.  For  simplicity,  we  consider  cold  ions  =  0).  Also,  we 

consider  a  weakly  collisional  plasma  in  the  sense  that  v  . ,  v  <<  ft  and 

ei  en  e 

v.  ,  v.  <<  ft.  where  v  .  represents  the  collision  frequency  between  species 

X  6  in  1  Hu 

a  and  b. 

The  basic  equations  used  in  the  analysis  are  continuity,  electron  and 
ion  momentum  transfer  and  current  conservation: 


dv. 

-l 

dt 


B)  -  v.  V. 
in-i 


m.  n 
1 


(3) 


V  •  J  -  V  •  [en(v.  -  vj]  -  0 


(4) 


where 


m  n  v  .  ( v 
e  e  eiv~e 


and  a  =  e,  i  for  electrons  and  ions,  respectively.  The  rest  of  the  symbols 
have  their  usual  meanings.  The  electron  fluid  is  assumed  to  have  an 
equilibrium  drift  velocity,  in  the  z-direction  and  is  non-uniform  in  the  x- 

A 

direction,  i.e.,  V  =  VA(x)  e  .  The  zero-order  current  is  expressed  as 
— U  U  2 

J0(x)  =  -  nQe  VQ(x)  ez  (5) 


where  we  take 


VQ(x)  =  VQ  exp  (-  x2/L2}  -  VQ(l  -  x2/L2}  (6) 

Here  we  have  used  a  parabolic  representation  of  the  finiteness  of  the 

current  along  the  x-axis,  which  is  an  approximation  for  the  normal 

distribution  with  a  half-width  of  L,,. 

w 

The  standard  procedure  is  followed  in  carrying  out  the  linear 
stability  analysis.  The  plasma  quantities  are  split  into  equilibrium  and 
perturbation  parts,  f  =  f^  +  <5f  ,  and  the  perturbed  quantities  are 
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assumed  to  vary  as  Sf  (x)  -  6f  (x)  exp  fik  y  +  ik  z  -  iaitl.  The  perturbed 

a  a  v  y  z 

ion  and  electron  equations  of  continuity  and  momentum  transfer  are  combined 
to  yield 


and 


5n. 


c2s{k2v  -  32/3x2) 

2  ~2 
“1  "  °i 


(7) 


veu2r1 

~2~2^  * 
k  v 
z  e 


r* 


(8) 


2  2 

where,  w,  =  to  +  iv.  ,  c  =  T  /m. ,  v  =*  T  /m  ,  v  =  \>  +  v  .  ,  ^  =  e6<f>/T 

1  in  s  ei  e  ee  e  en  ei  v  e 

and  to,,  =  to  -  k  V„(x) . 

2  z  0 

In  deriving  (7)  and  (8)  we  have  made  several  simplifying 

assumptions.  The  ion  temperature  is  neglected  (T^  =  o),  parallel  ion 

2  2  2  2 

motion  is  neglected  (w  >>  k  c  ),  k  /k  >>  v  /Q  is  assumed,  and  the 

ZS  z  6  6 

electrostatic  assumption  is  used  (SE  =  -  V6(j>).  We  further  make  use  of  the 
quasi-neutrality  assumption  (sn^^  *  5ne)  to  derive  the  nonlocal  mode 
structure  equation.  From  (7)  and  (8)  we  find  that 


♦  C(xH  -  0  (9) 

dx 


where 


Q(x)  «  [“  ky  +  W-  (to?  -  l)r], 
“l 


(10) 
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r  is  defined  in  (8),  and  we  have  written  the  variables  in  the  following 


dimensionless  form:  to  =  io/ft.t  x  =  x/p  ,  to,  =  <o  +  iv  , 

1  si  m 


"  “  *  k.,V 


va  =  W  ps  =  W  ky  =  kyps*  kz  =  kzps’  d  =  Lw/ps’  V0  =  W  and 
a  a  **  ^2  2 

VQ(x)  =  VQ(l  -  x  /d  ).  For  convenience,  we  will  drop  the  caret  over  the 


symbols . 


III.  ANALYTICAL  AND  NUMERICAL  RESULTS 


We  solve  (9)  numerically  for  the  parameters  appropriate  for  the 

auroral  ionosphere  and  obtain  the  nonlocal  growth  rate  of  the  collisional 

ion-cyclotron  modes  modified  by  finite  current  channel  width  effects. 

Before  presenting  these  results,  we  first  present  approximate  analytic 

2 

solutions  of  (9).  For  '<  »  \>  ai„,  we  write  (9) 

z  e  2 


Mf  +  [b  -  C  X2]^  =  0 


where 


I*  *y  <■?  -  Ul'  *  1  -^11 

Z 


o  v  Vn 

2  ,  i  e  0 


U:  -  i) 

0),  1  ' 


Equation  (12)  is  of  the  form  of  Weber's  equation  which  has  solutions 
determined  in  terms  of  Hermite's  functions.  The  eigenvalue  is  determined 


=  (2m  +  1  )C 
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Physically,  the  instability  arises  due  to  the  Doppler  effect  caused  by 
parallel  electron  streaming.  This  results  in  wave  growth  via  electron 
dissipation  (v>  )  when  the  electron  drift  speed  exceeds  the  parallel  wave 
phase  velocity  [Chaturvedi  and  Kaw,  1975]. 


The  nonlocal  growth  rate  is  given  by  [from  (15)] 


i  r  2  .  -i  veV0  1 1  / 2 


where  we  have  assumed  v^n  =  0  for  simplicity.  The  condition  for  marginal 
stability  (y  =  0)  is  therefore,  approximately, 


L  k 

d  =  —  =»  r  z 

ps  /2  V„k  ve^r 


if  2  -)  ve  ^0  ] 

2  iUp-1)—  1 


For  Vn  -  30,  v  =10  ,  k  -  0.5,  v.  -  0,  one  finds  L  *  Up  . 

0  e  y  m  w  Ks 

We  note  that  a  similar  stabilization  criterion  was  also  obtained  by 
Bakshi  et  al.  [1983]  for  the  collisionless  ion  cyclotron  instability.  The 

physical  interpretation  of  the  stabilization  of  the  mode  due  to  the 

finiteness  of  current  channel  width  is  as  follows.  The  finite  current 
channel  profile  considered  here  has  a  maximum  value  of  VQ  and  tends  to 

Vq  -*■  0  at  |X|  -»  d.  The  nonlocal  growth  rate  is  the  mode  growth  due  to  all 

the  regions  of  Vq(x)  that  are  sampled  by  the  wave  packet,  and  is  zero  for 
j  x |  >  d.  Clearly,  the  growth  rate  obtained  in  the  finite  width  channel 

case  is  reduced  from  the  case  in  which  the  current  sheet  would  be  infinite 
at  its  peak  value,  VQ.  Thus,  the  wave  packet  "sees"  an  effectively  reduced 
electron  drift  speed  in  the  finite-width  current  channel  case,  and,  for 
sufficiently  narrow  channels,  this  'effective'  drift  speed  may  not  be  large 


enough  to  exceed  the  parallel  wave  phase  speed  so  that  the  mode  is 


stabilized. 


The  nonlocal  wave  equation  (14)  was  solved  numerically  for  parameters 
appropriate  to  the  auroral  ionosphere.  The  growth  rate  (Y)  was  computed  as 


a  function  of  the  half-width  of  the  current-channel  (Lw).  Figures  2-5  show 

the  behavior  of  the  nonlocal  growth  rate  of  the  current-dr iven  collisional 

ion-cyclotron  instability  for  various  parametric  dependences.  For  all  the 

cases  considered,  the  transverse  wavelength  was  chosen  such  that  k  p  “0.5 

7  s 

and  the  parallel  wavelength  was  chosen  corresponding  to  the  maximum  growth 

-3 

rate.  In  Fig.  2  we  plot  Y/fl  versus  L^/p^  for  vg/fte  *  10  and 

v.  =0,  and  two  drift  velocities,  Vn/ca  =  20  and  30.  We  find  that  for 
m  ’os 

larger  drift  velocities,  complete  stabilization  occurs  for  narrower 
channels,  or  in  other  words,  the  instability  persists  for  smaller  width- 
channels  than  in  the  case  of  smaller  drift  velocity.  Figure  3  depicts  the 
dependence  of  the  Y/ih  versus  Lw/pg  as  a  function  of  electron  collision 
frequency  fv  1  for  Vn/ ca  =  30  and  v.  =  0.0.  Since  the  instability  is 
resistive  in  nature,  the  growth  rates  are  higher  for  larger  collisional 
frequencies.  Therefore,  for  mode  stabilization,  the  higher  the  electron 

collision  frequency  the  narrower  the  channel-width  required.  The  effect  of 
ion  collisional-damping  is  illustrated  in  the  Fig.  4  where  we  plot  '1/3. 

L 

-3 

versus  L  /p  for  Vn/c0  =  30  and  v  /fi  *  10  and  several  values  of  v.  /w.  . 
w  s  o  s  e  e  in  i 

The  inclusion  of  the  ion-damping  introduces  a  threshold  value  for  the 

electron  drift  velocity  for  the  instability  which  is  larger  than  the 

parallel  phase  velocity  of  the  mode.  Thus,  in  the  presence  of  ion 

collisions,  the  mode  can  become  stabilized  for  wider  current-channels  than 

in  the  case  when  v.  =*  0.  Finallv,  in  Fig.  5  we  plot  Y/Q.  versus  L  /p 
m  i  w  s 

for  different  mode  numbers  for  the  parameters  =  0,  ve/ng  ”  10  Vq/cs 
=  30.  The  lowest  order  mode  (m  =  0)  has  the  largest  growth  rate  and  higher 
modes  have  decreasing  growth  rates.  Thus,  the  higher  order  modes  are 
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stabilized  for  wider  current-channels  compared  to  the  m  -  0  mode.  We  note 
that  the  stabilization  width  is  Lu/p„  *  ^  in  Fig.  5,  which  is  in  agreement 
with  the  current-channel-width  computed  for  stabilization  from  the  analytic 
expression  (21). 

IV.  DISCUSSION 

We  find  that  the  finite  transverse  width  of  a  current-channel  has  an 

overall  effect  of  reducing  the  growth  rate  of  the  collisional  current 

driven  ion  cyclotron  instability  from  the  value  obtained  in  the  local 

approximation,  and  for  sufficiently  narrow  channels  can  stabilize  the 

instability.  For  typical  auroral  ionosphere  parameters,  i.e.,  Vq/cs  =* 

30,  v  /£!  -  10  ,  v./ft.  -  .02,  it  seems  that  the  complete  stabilization  of 

e  e  l  x 

the  instability  can  occur  for  current-channel  widths  L  /p  <  25  (or 
Ly  <  250  m  since  pg  -  10  m).  Observations  indicate  that  the  currents  often 
flow  in  sheets  with  thicknesses  in  ionosphere  on  the  order  of  -  km  [Burke 
et  al.,  1982].  Thus,  it  is  possible  that  this  instability  may  still  be 
operative  in  ionospheric  situations  provided  the  currents  are  intense 
enough  (so  that  the  electron  drift  velocity  is  above  the  threshold 
level).  These  results  are  qualitatively  similar  to  those  of  Bakshi  et  al. 
[1983]  who  considered  the  case  of  current  driven  collisionless  ion- 


cyclotron  instability. 


There  are  some  recent  laboratory  experiments  on  the  excitation  of  ion- 
cyclotron  instability  by  a  field-aligned  current  [Cartier  et  al.,  1985]. 
3y  varying  parameters,  one  may  be  able  to  scan  the  domain  of  current  driven 
EIC  instability  in  the  collisionless  and  collisional  domain  in  these 
experiments  and  thus  check  the  results  of  the  present  work  (collisional 
domain)  vis-a-vis  the  work  of  Bakshi  et  al.  [1983]  (collisionless 
domain).  We  note  that  the  theoretical  treatment  of  both  these  domains  for 
the  current-driven  EIC  instability  has  been  recently  carried  out  in  the 
local  approximation  by  Satyanarayana  et  al.  [1985]. 
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B. 

WHITE 

01  CY 

ATTN 

R. 

STAGAT 

MISSION  RESEARCH  CORP  . 

1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE,  NM  87106 
0 1 C Y  R.  STELLINGWERF 
0 1 C Y  M,  ALME 
01 CY  L.  WRIGHT 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  V A  22101 

0 1 C Y  ATTN  W.  HALL 
0 1 C Y  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1  23*10  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01 CY  ATTN  E.C.  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

0 1 CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

01 CY  ATTN  IRVING  L.  KOFSKY 


PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  4  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
0 1 CY  ATTN  FORREST  GILMORE 
01 CY  ATTN  WILLIAM  B.  WRIGHT,  JR 
0 1 CY  ATTN  WILLIAM  J.  KARZAS 
01 CY  ATTN  H.  ORY 
0 1 CY  ATTN  C.  MACDONALD 

RAND  CORPORATION,  THE 
15450  COHASSET  STREET 
VAN  NUYS ,  CA  91406 

01 CY  ATTN  CULLEN  CRAIN 
0 1 C Y  ATTN  ED  BEDROZI AN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01 CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

0 1 CY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

01 CY  ATTN  LEWIS  M.  LINSON 
0 1 CY  ATTN  DANIEL  A.  HAMLIN 
0 1 CY  ATTN  E.  FRIEMAN 
0 1 C Y  ATTN  E.A.  STRAKER 
0 1 C Y  ATTN  CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
0 1 C Y  J.  COCKAYNE 
0 1 C Y  E.  HYMAN 


PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  3027 
BELLEVUE,  WA  98009 

0 1 C Y  ATTN  E.J.  FREMOUW 
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SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 
0 1  C  Y  ATTN  J.  CASPER 
0  1  C  Y  ATTN  DONALD  NEILSON 
01CY  ATTN  ALAN  BURNS 
0 1 C Y  ATTN  G.  SMITH 
0 1 C Y  ATTN  R.  TSUNODA 
01 CY  ATTN  DAVID  A.  JOHNSON 
0 1 C Y  ATTN  WALTER  G.  CHESNUT 
0 1 C Y  ATTN  CHARLES  L.  RINO 
0 1 C Y  ATTN  WALTER  JAYE 
0 1 C Y  ATTN  J.  VICKREY 
01 CY  ATTN  RAY  L.  LEADABRAND 
0 1 C Y  ATTN  G.  CARPENTER 
0 1 C Y  ATTN  G.  PRICE 
0 1 C Y  ATTN  R.  LIVINGSTON 
01CY  ATTN  V.  GONZALES 
0 1 C Y  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01 CY  ATTN  W.P.  BOQUIST 

TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
0 1 C Y  ATTN  R.  K.  PLEBUCH 
0 1 C Y  ATTN  S.  ALTSCHULER 
01 CY  ATTN  D.  DEE 
0 1 C  Y  ATTN  D/  STOCKWELL 
SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MA  01803 
0 1 C Y  ATTN  W.  REIDY 
0 1  C Y  ATTN  J.  CARPENTER 
0  1  C Y  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

0 1 CY  ATTN:  N,  ZABUSKY 

DIRECTOR  OF  RESEARCH 
U.S.  NAVAL  ACADEMY 
ANNAPOLIS,  MD  21402 
02CY 

CODE  1220 
01CY 
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